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We calculate the hydrogen molecule ion from the two particle Schròdinger equation. Therefore 
a very simple two particle basis set is chosen. We suggest this ansatz to be used to solve the "two 
electron one phonon" three particle wave-function of a BCS superconductor. Possibly it can give 
hints for high temperature superconductors. 



I. INTRODUCTION 

Ï/V 

Quantum mechanical calculations of rotational, vibrational and electrònic excitations in molècules are usually done 
in the Born Oppenheimer approximation [1]. The electron-phonon system of crystals is usually calculated in this 
approximation, too. Coupling of vibrations and electrònic excitations is calculated due to perturbation theory. In 
this article we take the hydrogen molecule ion as a very simple example of a system with rotational, vibrational 
) I ■ and electrònic excitations and suggest a way of calculating it without using the Born Oppenheimer approximation. 
Therefore we will use two particle wave-functions. The same concept may be used to calculate electron-phonon 
systems of crystals. 
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II. CÀLCUL ATION 



We use the Schròdinger equation 



O '■ H\<p>=e\tp> (1) 

in 

for the hydrogen molecule ion. The problem has three partides, but it can be reduced to a two particle problem in 
a similar way, as the hydrogen atom can be reduced from a two particle problem to a one particle problem [1]. The 
Hamiltonian has the form 

S : ff = -^A a -^A b -^A c + -^ + -^ + -^. (2) 

2m Q 2m b 2m c \r a - n\ \n - r c \ \r a - r c \ 

3 ! 

A a , Af, and A c are the Laplace operators with respect to the positions of the partides r a , r b and r c , respectively. 

• • . We introduce the center of mass 
> ■ 

'/^ 1 _ m a r2 + m b rj + m c rt 

\2 ■ m a + m b + m c 

and get new variables 



rí = n + Tl,n = rï + ^. , (4) 

_> _ m a ri + m b Tb _ m a (rí + ~È) + m b ij2 + ^) /_>. 

7 ' c — — ~Z ■ v*) 

m c m c 

We are only interested in solutions, which do not depend on the motion of the center of mass, therefore we use 
Tt = 

r a = n , r b = r 2 , r c = . (6) 

m c 
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Let m a — m e be the electron mass, mi, = m c = m p be the proton mass, e a = — e the electron charge and et = e c = e 
the proton charge. With m a <C m c and mj, = m c we get 

rï = -r$. (7) 

This is not a Born Oppenheimer approximation, it is just a simplification of the calculation. The value of r c can 
be wrong by 0.0005 of usual distances in the hydrogen molcculc ion. Within our calculation of the matrix elements 
this will be a small error with respect to the integration errors from Monte-Carlo integration. 

Using the new variables r? and fi we get a two particle Hamiltonian 

__ h K e 2 e 2 e 2 

= ~2^i x ~^ v 2 + 2M _ if?TWj ~ ' ( ) 

Ai, A2 are the Laplace operators with respect to ff and fi, respectively. Due to constrains the mass moved by 
changing fi is 2m p . With this Hamiltonian we calculate the matrix elements for the Schrodingcr equation for a chosen 
two particle basis-function set. We use the linear combination ansatz and the variation theorem 

<fi(ri, ri) = ^íCí<Pí(TÍ, fi) . (9) 
From the expectation value e of the energy 

< ip\H\ip >= e < ip\<p > , (10) 

we get 

J dV 1 dV 2 T, ij c j ip* j (ft,fÍ)Hc i <fi i (ft,fÍ) = J dV 1 dV 2 T, ij c j <fi*(ft,fÍ)ec i <fi i (ft,fÍ) , (11) 

and write it as 
with 

Ha = J dV 1 dV 2 tp l Htp 3 , S t3 = j dVxdVïViipj . (13) 

The integrals must be calculated over two 3D-volumes, as ipi are two particle basis-functions. From the variation 
theorem the energy e is a minimum with respect to all parameters Cj 

dej YlijCjSijCi 

With this ansatz the Schrodingcr equation looks likc 

^iHijCi = e YiiSijCi . (15) 

This is a general eigenvalue problem which can be solved numerically. Also the integrals Hij and SV, can be solved 
numerically, for example by Monte Cario integration. 

III. BASIS-FUNCTIONS FOR SOLVING THE TWO PARTICLE SCHRÒDINGER EQUATION 

The choice of the two particle basis-functions <Pi(fí,fÍ) will bc cvaluatcd now. If two particle electron wave- 
functions are needed, often the ansatz (p pq (fï,fi) = ' l P' p {fï)4>'q{^ï) is used. In our case this ansatz will be very bad, 
because one particle is a (quasi) proton, the other an electron. The motion of the electron is not independent from 
the motion of the protons, as the electrons will have highest probability density near the protons. A simple approach 
could be 



Ppq(rí, fi) = k pq 1pp{fí) (<f>q{fÍ ~ fí) + S<f> q (fÍ + fi)) . 



(16) 
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Eigenvalue [a.u.] 


coo 


ClO 


C20 


C30 


-0.521112 


-0.99955 





-0.0299897 





-0.520852 





-0.999777 





-0.0211105 


-0.505818 


-0.198128 





0.980176 





-0.505597 





0.190136 





-0.981758 



Table I: Results of sample calculation 



We will use ip p and <fi q as radial symmetric functions. The combination of the indices p and q correspond to the index 
i — > pq in (15). k pq is the normalization constant and will be calculated from the overlap integrals J dVidV2(p P q(p pq = 1. 
(^9(^2 — rt) + s(f> q (r2 + rí)) describes the electron part of the basis-functions. The highest probability density of the 
electron is near both protons, s is +1 or — 1 characterizing the bonding and anti-bonding wave-functions of the 
electron [1]. We will only use s = 1 for bonding wave-functions and use simple basis-functions for 4> q and ij) p 

M?) = f p {r)e- w ^{\t\-rot ) 

<fo{~r*) = e- w " l "^"" ( ^^ 2 . (17) 
r is the half of the proton-proton distance. As functions f p we use: 

fo(T) = l , f 1 (^) = -^- , f 2 (j>) = (\r\-r ) , h(?) = f2(?)fi(?) , 



with T = y y J . (18) 

ipo should be similar to the wave-function for l = , n = (rotational and vibrational ground stage), tpi is similar 
to l = 1 , n = 0, tp2 is similar to i = , n = 1, ^3 is similar to Z = 1 , n = 1. Vo is constructed as ground wave-function 
of a harmònic oscillator, -02 is the first excitation of a harmònic oscillator with respect to the radial dependency (ro 
is the mean distance of the protons from the center of mass) [1]. describes the spheric harmonies for 1 = 1. 

4>o(~^) is a Gaussian type ground wave-function for the electron. With this basis-functions we choose approximate 
parameters in atòmic units 

VJproton — 18 > 

Welectron = 0.45 , (19) 
r = 1.0. 

This parameters can be calculated by searching the minimum of the ground state energy from equation (15) with 
respect to w pro ton, Welectron and tq. This is not done within this article, as the computational power was not available. 
We only did some test calculations to check, if we are near the minimum. The parameters are taken from known 
properties of the hydrogcn moleculc ion [2] . 

The results for the four cigcnvalues and eigenvectors are shown in table I. They are not very good compared to 
other calculations, but of cause the choice of the basis set is very simple. The ground state wave-function is close to 
</?oo, the first excitation is the rotational excitation with a wave-function close to </?io- The second excitation is the 
vibrational excitation with a wave-function close to </?20- The calculated energy for the first rotational excitation is 
0.00026 a.u. = 0.007 eV. The calculated energy for the first vibrational excitation is 0.0153 a.u. = 0.41 eV. The vàlues 
can be compared to literature [2] . The energy of the vibrational excitation is calculated to high. The literature value 
is about 0.27 eV. One can see from table I, that the eigenfunctions arc not very good. In the linear combination also 
c o diffcrs from zero. 



IV. DISCUSSION 



The suggested basis-functions (16) can be used to calculate the two particle problem of the hydrogen molecule 
ion. As result we get two particle wave-functions for the ground state, the first rotational excitation and the first 
vibrational excitation. 
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The chosen basis-functions are very simple functions and the basis set is very small. They can of cause be improved 
in a number of ways. Espccially for the clcctron part of the wave-function we did not even include basis-functions for 
electrònic excitations. 

Additionally the separation (equation: 16) does not allow electron phonon interaction, as the electron wave-function 
does not depend on the protons. One will need electron basis-functions, which include the proton-proton distance 2|f2 1 
to include electron-phonon interaction. This is not a principle problem, as the separation is only used to construct 
simple basis-functions. The separation is not needed for the calculation in any way. With enough computational 
power it should be possible to use a subset of a complete set of basis-functions. With the use of such a subset, no 
artificial approximations are needed to solve the Schròdingcr equation and the precision can be increase by increasing 
the number of basis-functions in the subset. 



V. OUTLOOK 



This ansatz can be transferred to crystal lattices. The aim would be to calculate the electron-phonon system in 
a one electron, one phonon (more general: vibrational excitation) approximation ab initio. In a next step a "two 
electron one phonon" approximation could be used to calculate the electron-phonon system of a BCS superconductor. 
With the choice of a subset of a complete set of three particle basis-functions no additional approximations will be 
needed (but the limited number of basis-functions in the subset). It would be quite interesting to see, how this 
wave-functions look like. 

Possibly this ansatz can give new ideas, how to describe high temperature superconductors. 

VI. RESOURCES 

The source code of this calculation is available [3]. Please contact the author in case you want to work on the tòpic. 
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